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ABSTRACT 

We confirm the UHEGR horizon estabhshed by the Pierre Auger Observatory using the heteroge- 
neous Veron-Getty Veron (VCV) catalog of AGNs, by performing a redshift-angle-IR luminosity scan 
using PSGz galaxies having infrared luminosity greater than 10^° Lq. The strongest correlation - 
for z < 0.016, ij} = 2.1°, and Lir > 10^'^'^ Lq - arises in fewer than 0.3% of scans with isotropic 
source directions. When we apply a penalty for using the UHEGR energy threshold that was tuned 
to maximize the correlation with VGV, the significance degrades to 1.1%. Since the PSGz catalog is 
complete and volume-limited for these parameters, this suggests that the UHEGR horizon discovered 
by the Pierre Auger Observatory is not an artifact of the incompleteness and other idiosyncrasies 
of the VGV catalog. The strength of the correlation between UHEGRs and the nearby highest-IR- 
luminosity PSGz galaxies is stronger than in about 90% percent of trials with scrambled luminosity 
assignments for the PSGz galaxies. If confirmed by future data, this result would indicate that the 
sources of UHEGRs are more strongly associated with luminous IR galaxies than with ordinary, lower 
IR luminosity galaxies. 

Subject headings: cosmic rays, AGN, x-rays, GRB 



1. INTRODUCTION 

The Pierre Auger Observatory has reported 
(|The Pierre Auger Collaboration! [20071 [20081 ) a cor- 
relation between the arrival directions of the highest 
energy cosmic rays, and the positions of galaxies in 
the Veron-Cetty Veron "Catalog of Quasars and Active 
Gala ctic Nuclei" (12th Ed.) (jVeron-Cettv fc VeronI 
|2006( ) (VCV below). Of the 27 cosmic rays above 57 
EeV recorded by the Pierre Auger Observatory prior 
to Aug. 31, 2007, twenty are within 3.2° of a VCV 
galaxy with z < 0.018 (about 75 Mpc). Restricting 
to \b\ > 10° to avoid regions in which galaxy catalogs 
are incomplete due to Galactic extinction, leaves 22 
UHEGRs of whic h 19 are within 3.2° of a VCV galaxy 
(jZaw et al.l I2O1O0 . The existence of a correlation be- 
tween UHEGRs and VGV galaxies clarifies the nature 
of ultrahigh energy cosmic rays by establishing that 
UHEGRs are extragalactic in origin and that they have 
an energy- dependent horizon consistent with the GZK 
prediction. However the VCV catalog is not complete, 
homogen eous, or pure, and its comple t eness drops with 
redsh ift ([Veron-Cettv fc VeronI [200l iZaw et all [2001 
|2010[) . Thus VGV is not suitable for statistical analyses 
as would be required for establishing quantitative 
bounds on the contribution of various possible sources 
to the population of UHEGRs, and its suitability for 
establishing the GZK horizon can be a source of concern. 

Here, we study the correlation between UHEGRs and 
infrared luminous galaxies with two aims(3 First, we 
use the IRAS point source catalog with redshifts (PSGz, 
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iSaunders et al.|[2000 D. which is homogeneous and com- 
plete to beyond the GZK distance, to replicate the Auger 
result of a UHEGR horizon and correlation with ex- 
tragalactic matter. Second, we investigate whether the 
strength of a galaxy's correlation with UHEGRs is re- 
lated to its infrared luminosity, Ljr. In the absence of an 
active nucleus, infrared luminosity is a good proxy for 
star formation, and thus galaxies with strong IR emission 
should have an enhanced concentration of massive stars, 
whose lifetimes are short on the time scale of star forma- 
tion episodes. A number of proposed UHEGR sources 
are associated with the death of massive stars or birth 
of their progeny, e.g., collapsars, GRBs, or magnetars, 
so the probability of such events occurring in a given 
galaxy could increase roughly in proportion to iir of the 
galaxy. These considerations motivate a study of the 
correlation between UHEGRs and high Li^ galaxies. Of 
course, some high Lir galaxies host AGNs, so discovering 
that correlations are enhanced with a high Lji- sample 
does not by itself demonstrate a non-AGN source popu- 
lation and further investigation is needed. In particular, 
low-redshift GRBs occur primar ily in relatively rare low 
metallicity star-forming regions (|Stanek et al.l[2nM ). es- 
pecially found in dwarf galaxies, so a UHEGR correlation 
with high luminosity IR galaxies would not necessarily 
be indicative of GRB sources. The correlation between 
A GASA UHEGRs and PSG z galaxies was investigated 
bv'Smialk owski et al.l ()2002[ ) and the cluste ring of Auger 
UHEGRs was compared to that of PSGz bv lCuoco et al.1 
(2008.) . 

We use the technique employed by Auger, of scan- 
ning over the maximum angular separation, ip, and 
over the maximum redsh ift of the galaxies (also see 
iTinvakov fc Tkachevl[200lb . We also scan on the thresh- 
old infrared luminosity of galaxies. Since only events 
with E > 57 EeV have been published, we do not do a 
full scan on energy, but we do a partial scan in one di- 
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rection (toward higher energy) in order to estimate the 
penalty associated with adopting the 57 EeV threshold. 
To assess the significance of the correlation compared to 
isotropy, as relevant to establishing the existence of a 
UHECR horizon, we perform the same scans on mock 
isotropic UHECR datasets. To test the significance of a 
possible correlation with infrared luminosity, we compare 
to the results of scans using mock catalogs obtained by 
scrambling the luminosities of the PSCz galaxies. 

2. SCAN ANALYSIS WITH INFRARED LUMINOUS 
GALAXIES 

IRAS is an all-sky survey of infrared galaxies, and 
the PSC z catalog provides t he red shifts of IRAS point 
sources. iSanders fc Mirabel (|1996l ) give the relationship 
between the IRAS observed fluxes at 12, 25, 60 and 
lOO^m in W/m^ ifi2, etc.), and the total 8-1000/^m in- 
frared luminosity; 

Li, = AttDI 1.8x10-14[13.48/i2+5.16/25+2.58/6o+/ioo]^ 

(1) 

where is the luminosity distance. Figure [T] shows the 
distribution of log[Lii/Lo] in the entire PSCz catalog 
(solid), and, for later reference, for z < 0.016 (dotted). 
We implement the PSCz mask by a simple Galactic lat- 
itude cut \b\ > 10°, the recommended conservative cut 
to avoid poorly samp led regions near the Galactic plane 
(jSaunders et al.ll2000[) . 

We perform a scan on the angular separation, maxi- 
mum redshift and minimum infrared luminosity of PSCz 
galaxies and UHECRs with |&| > 10°, covering the ranges 
V; = 0.1° - 5°, z„,ax = - 0.024, and log[Lir/Lo]min = 
10 - 11, in steps of A?A = 0.1°, Az = 0.001, and 
Alog[Li,./L0]min = 0.1. For each combination of 
ip, and log[Lir/L0]min we find the number of correlated 
UHECRs, fccorr, and compute the associated probabil- 
ity measure P. P is the cumulative binomial proba- 
bility to find fccorr correlated events out of A^cR cosmic 
rays, given the fractional sky coverage p of PSCz galaxies 
within the given redshift and luminosity bounds, p is the 
exposure- weighted projection of disks of radius ip around 
each galaxy in the sample, normalized to the total Auger 
exposure; it must be recomputed for each bin in Zmax, 
V', and log[Lir/i0]min- The minimum value of P in the 
scan analysis is found for z,„ax = 0.016, log[Lir/i0]min = 
10.5, and ip = 2.1°. With these parameters, there are 
632 PSCz galaxies, p = 0.1611, and 13 out of 22 UHE- 
CRs correlate. The associated probability measure is 
P = 5.7 X 10~^. The probability of finding as good 
or better correlation by chance from an isotropic source 
distribution is (2.68 ± 0.07) x 10~^. This is measured 
by creating 5 x 10^ catalogs of 22 random UHECR direc- 
tions with \b\ > 10° weighted by the Auger exposure, and 
using these to perform the scans with the PSCz galaxy 
sample rather than real Auger data. 

The choice of energy cut {E > 57 EeV) made by Auger 
in releasing the data was tuned to maximize the corre- 
lation of UHECRs with VCV AGN. While this analysis 
is independent of the VCV analysis performed by Auger, 
VCV and PSCz galaxies both trace the same underly- 
ing matter density field. It is thus reasonable to worry 
that the choice of energy cut will bias our results. In 
order to estimate a penalty for this, we include a scan 
in energy threshold in one direction, toward higher en- 
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Fig. 1. — Luminosity distribution for the full PSCz catalog 
(solid) and within z < 0.016 (dotted). The vertical line at 10.5 in- 
dicates the luminosity threshold for maximally significant UHECR 
correlations. Note that there are very few luminous PSCz galaxies 
at low redshift. 
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Fig. 2. — Infrared luminosity vs. redshift for PSCz galax- 
ies. The vertical and horizontal lines at indicate the values of 
log[Lir/L0]min and ^max for which UHECR correlations are max- 
imum. 

ergy. We find that the minimum value of P occurs for 
the original energy cut of 57 EeV and the same values 
of Zmax J and log [Lii./L0] mill- We then use 5 x 10** 
catalogs of 22 random UHECR directions in new scans 
that include the energy scan in one direction. We find 
that the probability of finding as good or better corre- 
lation by chance from an isotropic source distribution is 
(1.10 ± 0.05) X 10~^. In other words, the significance of 
our result degrades from 0.3% to 1.1% when we include 
the energy penalty. 

As can be seen from Figure^ the PSCz catalog is com- 
plete to redshift 0.024 for log[Lii./Lo] > 10.6, to 0.02 
for log[Lir/L0] > 10.5, and for z < 0.016 it is com- 



3 



plete down to log[Lii/LQ] < 10.3. The scan minimum 
(log[_Lii./_L0]min = 10.5, Zmax = 0.016) is thus Comfort- 
ably inside the domain of PSCz completeness and is thus 
established with a volume-limited sample. Although the 
whole scan range includes regions for which the PSCz 
catalog is not complete, the simulations measuring the 
chance probability for isotropically distributed UHECRs 
to give as low or lower a P value have the same in- 
completeness, and therefore give a reliable estimate of 
the significance. Therefore, the analysis presented here 
suggests that the arrival directions of UHECRs corre- 
late with nearby extragalactic objects and establishes the 
UHECR horizon. 

3. INFRARED LUMINOSITY AND UHECR CORRELATION 

One next wants to know if this correlation of UHE- 
CRs with luminous PSCz galaxies is stronger than if the 
galaxies were selected from the PSCz catalog without re- 
gard to their IR luminosity. To find out, we follow the 
same scan procedure as above, but scrambling the val- 
ues of log[Lir/I/0] first. The P values obtained in 700 
such trials are shown in the upper panel of Figure [3] Six 
percent of the cases have a P value as low as in the real 
scan using the highest luminosity IR galaxies. 

Two of the VCV galaxies that correlate with pre- 
prescription UHECRs are also in the PSCz catalog. Thus 
the formulation of the Auger prescription might indi- 
rectly bias the present study of UHECR-LIRG correla- 
tions. We can remove any such bias in estimating the 
chance probability, at the price of obtaining only a bound 
on the significance, by removing those two UHECRs from 
the UHECR data set. Since that explicitly eliminates 
some correlations, we find an upper bound on the a pri- 
ori chance probability. The lower panel of Figure [3] shows 
the distribution of P values for the reduced UHECR data 
set and scrambled IR luminosity values; the vertical line 
is the P value found using the reduced UHECR data set 
and scanning on the true IR luminosity. The observed P 
value is lower than in the scrambled data sets in 11.7% of 
the cases, allowing us to bracket the chance probability 
for the UHECR-LIRG correlation not being Lji. depen- 
dent to be in the range 6-12%. This result suggests that 
UHECRs are more strongly correlated with highly lumi- 
nous PSCz galaxies than lower luminosity galaxies, but 
it is not statistically significant and remains to be con- 
firmed with future data. 

Table 1 gives the name. Right Ascension and declina- 
tion of the PSCz galaxies that correlate with UHECRs 
for the scan parameters that minimize P {z < 0.016, 
^ < 2.1° and logtLir/L©] > 10.5), followed by the year, 
Julian day, and energy of the correlated UHECR. This is 
followed by the angular separation between galaxy and 
UHECR, and the redshift and derived infrared luminos- 
ity of the galaxy. In some cases, PSCz reports mul- 
tiple values of Li,-; in such cases these are listed sep- 
arated by a comma. Thirteen distinct PSCz galaxies 
correlate with 13 different UHECRs within 2.1°, with 
two UHECRs correlating with more than one PSCz 
galaxy and two PSCz galaxies correlating with more 
than one UHECR: 2MASX J17544125-6054404 {IRAS 
17501-6054) and ESQ 270- G 007 {IRAS 13203-4317). Of 
course, a given cosmic ray can have only one source, and 
the presence of two source candidates within the scan- 
ning prescription is a useful reminder of the limitations 
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Fig. 3. — Upper panel: histogram of P values obtained in 700 
scans of UHECRs and PSCz galaxies, but for scrambled values of 
log[Lir/LQ]; the P value for the true log[Lir/-L0] is shown as a 
vertical line. Lower panel: same, but with 2 UHECRs correlating 
with VCV galaxies removed. 

of the scanning method. 

It must be emphasized that a correlation in the scan 
procedure does not imply a correct identification of the 
CR source: the true source for that cosmic ray may not 
even be in the catalog or may be beyond the Zmax which 
minimizes P. Only a significant excess of correlations 
compared to chance, for the ensemble of correlations, is 
evidence for a concentration of sources in the c atalog. 
The Generalized Maximum Likelihood method (jFarraii 
l2008f) . which allows the GZK redshift information to be 
imposed on an event-by-event basis, is distinctly supe- 
rior to the scan method in this regard (Busca and Far- 
rar, 2009, in preparation) but we have adopted the scan 
method for the present work, to allow direct comparison 
to the Auger analysis with the VCV catalog. Use of the 
new Generalized Maximum Likelihood method should 
significantly improve the sensitivity of future correlation 
studies (Zaw, Busca, Farrar, 2010, in preparation). Since 
PSCz is complete for the redshift and L,,- range of in- 
terest, and outside its mask the PSCz catalog does not 
suffer from the inhomogeneity issues plaguing the VCV 
catalog, detailed statistical studies to gain further insight 
into the origin of the correlation will be possible. 

Both AGN activity and star formation lead to a high 
Lii- after reprocessing by dust, and indeed 5 of the 13 cor- 
related PSCz galaxies have opt ically ident ified active nu- 
clei by the criteria described in lZaw et al.l ((2009). This is 
a much larger portion than for randomly chosen high IR 
luminosity galaxies, since studies show that 4% of galax- 
ies with log[Lir/L0] in the r ange 10 — 11 are Seyferts , 
i.e. confirmed optical AGNs (jSanders fc MirabellflQOl . 
Only ab out half of a ll active nuclei can be identified in the 
optical (|Reviglio fc Hclfand 2006), so an important ques- 
tion is whether any of the 7 correlated PSCz galaxies that 
do not appear in VCV have active nuclei as well. Upcom- 
ing X-ray observations with Chandra of these galaxies 
will help settle this question. For the moment, we can- 
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TABLE 1 

Properties of correlating PSCz galaxies 



IR Galaxy RA Dec Auger Yr, Day E r z Li, AGN? VCV Name 





(J2000) 


(J2000) 




(EeV) 


(deg) 










NGC 1346 


52.55312 


-5.54282 


2006, 296 


69 


1.04 


0.013 


10.52 


SI 


SDSS J03302-0532 


NGC 1204 


46.16654 


-12.34082 


2007, 145 


78 


1.56 


0.014 


10.91 


H2 


NGC 1204 


ESO 565-G006 


141.80370 


-19.53184 


2007, 84 


64 


1.83 


0.016 


10.81 






NGC 424 


17.86485 


-38.08423 


2005, 54 


84 


0.37 


0.012 


10.62 


SI. 5 


NGC 424 


NGC 7130 


327.08014 


-34.95164 


2007, 193 


90 


2.00 


0.016 


11.35 


S2 


IC 5135 


NGC 4945 


196.36290 


-49.46769 


2005, 81 


58 


2.02 


0.002 


10.79,11.03 


S2 


NGC 4945 


ESO 270-G007 


200.81441 


-43.54422 


2006, 299 


69 


1.75 


0.013 


10.55 












2007, 69 


70 


0.49 










NGC 5128 (Ceil A) 


201.36003 


-43.01828 


2007, 69 


70 


0.91 


0.002 


10.52.10.62 


BLLac 


NGC 5128 


IC 4523 


226.29191 


-43.50967 


2004, 343 


63 


1.46 


0.016 


10.57 






2MASX J17544125-605440 


268.67473 


-60.91092 


2004, 339 


83 


0.14 


0.016 


10.52 












2006, 55 


59 


0.54 










NGC 7648 


350.97771 


9.66911 


2006, 185 


83 


0.93 


0.012 


10.55 






IC 5179 


334.03740 


-36.84463 


2005, 295 


57 


1.68 


0.011 


11.16 






IC 5186 


334.69389 


-36.80152 


2005, 295 


57 


2.05 


0.016 


10.65 







Note. — This table lists key properties of the IRAS PSCz galaxies with \b\ > 10° that arc within 2.1° of an Auger UHECR. The arrival year 
and Julian day, and the energy of the correlated UHECR are given for each galaxy, followed by the separation between the IR galaxy and UHECR 
in degrees, and the redshift and derived infrared luminosity of the galaxy; when multiple values of Ljr are given in PSCz, all are listed, separated 
by a comma. For galaxies appearing in the VCV catalog, the ZFG classification and VCV name for the galaxy is also given. 



not determine whether the UHECR correlation with IR 
luminous galaxies may result from AGNs they contain, 
or may represent evidence for a second class of sources. 



4. CONCLUSIONS 

The scan analysis performed here shows a signifi- 
cant correlation between the published Auger cosmic 
rays and nearby moderately-to-highly luminous infrared 
galaxies. This provides a valuable independent corrob- 
oration of the UHECR horizon with a homogeneous, 
volume-limited catalog. The correlation is maximized 
for z < 0.016, ^ = 2.1° and log[Lir/i0]min = 10.5. The 
probability that isotropic sources would show such a high 
degree of correlation by chance is 0.3%. The PSCz cata- 
log is complete to 2; = 0.024 for logfLir/L©] > 10.5, thus 
we establish the UHECR horizon with a volume-limited 
catalog, removing the concern that the horizon found in 
the Auger analysis with VCV galaxies could be an ar- 
tifact of VCV's i ncreasing incompleteness with distance 
(|Zaw et al.ll20T0l) . 

The angular separation maximizing the correlation in 
the PSCz luminosity scan, ip = 2.1°, is lower than the 
ij; = 3.2° optimizing the VCV scan. The relationship be- 
tween ip and the spectrum of magnetic deflections is sub- 
tle and depends on the source sample size and complete- 
ness, since ip is determined by minimizing P. Nonethe- 
less, this difference may be evidence that the incomplete- 
ness and inhomogeneity of VCV exaggerates the angular 
separation between source and UHECR arrival direction, 
and thus one must be cautious about inferring the typical 
magnetic deflection at this time. 

The correlation observed between the published Auger 



UHECRs and highly luminous PSCz galaxies is stronger 
than between UHECRs and randomly chosen PSCz 
galaxies within the same redshift and angular region, in 
6-12% of the trials. If confirmed by future data, such 
a correlation between UHECRs and high-IR-luminosity 
galaxies would show that luminous IR galaxies have an 
increased abundance of sources of UHECRs compared to 
randomly chosen galaxies in the PSCz catalog. It does 
not, however, resolve the question of whether UHECRs 
are produced by AGNs or some other mechanism, per- 
haps associated with the demise of stars as in GRBs or 
magnetars, since many of the correlated infrared galax- 
ies are known to contain AGNs. A program of Chandra 
observations is underway to find out which of the other 
correlated high-IR-luminosity galaxies may also have ac- 
tive nuclei. 

We note that new Auger data (since the published 
2007 Science results) show a weaker correlation with the 
VCV AGN (Hand 2010), highlighting that the published 
UHECR dataset may be too small to place robust con- 
straints on cross-correlations with extragalactic sources. 
It is thus especially important to repeat the present anal- 
ysis with a new and larger dataset. 
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berg for helpful input, as well as the anonymous ref- 
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